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ABSTRACT 

We detect the weak gravitational lensing distortion of 450,000 background galaxies 
(20 < R < 23) by 790 foreground galaxies (R < 18) selected from the Las Campanas 
Redshift Survey(LCRS). This is the first detection of weak lensing by field galaxies of 
known redshift, and as such permits us to reconstruct the shear profile of the typical field 
galaxy halo in absolute physical units (modulo H ), and to investigate the dependence 
of halo mass upon galaxy luminosity. This is also the first galaxy-galaxy lensing study 
for which the calibration errors due to uncertainty in the background galaxy redshift 
distribution and the seeing correction are negligible. Within a projected radius of 200 
h^ 1 kpc, the shear profile is consistent with an isothermal profile with circular velocity 
v c = 164 db 20 km s _1 for an galaxy, consistent with the typical circular velocity for 
the disks of spirals at this luminosity. This halo mass normalization, combined with the 
halo profile derived by Fischer et al. (2000) from a galaxy-galaxy lensing analysis of the 
Sloan Digital Sky Survey, places a lower limit of (2.7 ± 0.6) x 1O 12 /i _1 M0 on the mass 
of an L* galaxy halo, in good agreement with the satellite galaxy studies of Zaritsky et 
al. (1997) Given the known luminosity function of LCRS galaxies, and assumption that 
M cx L@ for galaxies, we determine that the mass within 260/i~ 1 kpc of normal galaxies 
contributes O, = 0.16±0.03 to the density of the Universe (for (3 = 1) or £1 = 0.24±0.06 
for (3 = 0.5. These lensing data suggest that 0.6 < (3 < 2.4 (95% CL), only marginally 
in agreement with the usual (3 ~ 0.5 Faber-Jackson or Tully-Fisher scaling. This is the 
most complete direct inventory of the matter content of the Universe to date. 
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1. Introduction 

Gravitational lensing provides the only non-dynamical measure of the mass distribution of a 
galaxy halo, and it thus enables a crucial means of quantifying the dark matter halos that are now 
an essential component of all theories of cosmology and galaxy formation. A reliable measure of the 
masses of galaxies and their halos also provides a significant lower limit to the density parameter 
S7. Spiral galaxy halos have traditionally been probed via rotation of stars and gas, but these 
tracers are not available outside radii of ~ 30 kpc, though their velocities are consistent with with 
an isothermal halo extending beyond this distance. At larger radii, dynamical test particles are 
scarce: Zaritsky et al. (1997) use satellite galaxies to trace the halo mass profile to « 200 kpc. 
Because there are only « 1.2 satellites per spiral, the masses of individual galaxies are poorly 
constrained, but by "stacking" the signals from an ensemble of galaxies, an accurate measure of 
the mean halo is obtained. A further complication of such dynamical studies of the outer halo is 
that the test particles are not virialized, and determination of the halo mass requires some model 
of infall. With such modelling, Zaritsky et al. (1997) estimate a mass of ~ 2 x 10 12 Mq within 
200 kpc of an L* spiral galaxy. 

The gravitational lensing approach uses photons from background galaxies (bggs) as the test 
particles, enabling halo measurements to large radii. The greatest problem with lensing is the 
extremely weak shear (< 3%) associated with the mass of a single galaxy. This is more than a 
factor of 10 below the noise due to intrinsic shape variations in a single background galaxy. As 
with the satellite galaxy study, the weak lensing approach must therefore overcome the poor S/N 
per galaxy by using large numbers of lens-source pairs, measuring the average galaxy halo. The 
weak lensing measurement, however, requires no modeling of the dynamical state and can give a 
non-parametric estimate of the halo profile. The first attempted detection was Tyson et al. (1984), 
and over the last five years several groups have verified the existence of galaxy-galaxy lensing based 
upon analysis of a small number of deep fields (Brainerd, Blandford, & Smail 1996; Dell' Antonio 
&; Tyson 1996; Griffiths, Casertano, & Ratnatunga 1996; Hudson et al. 1998). Interpretation of 
these results is problematic not only because of poor S/N, but also because the redshifts of the lens 
galaxies (and hence their distances and luminosities) are not known individually. Interpretation 
relies upon positing a model for the joint distance/luminosity distribution of the lenses and sources, 
and a model for dependence of halo mass and size upon luminosity, then projecting these models 
down to a function of angular variables to be compared with the observations. In such deep fields, 
even the statistical distributions of lens and source distances are currently ill-determined. 

More recently Fischer et al. (2000, F00) have measured a galaxy-galaxy lensing signal at high 
significance from preliminary Sloan Digital Sky Survey data. These data are relatively shallow 
(foreground galaxies 16 < r' < 18, background galaxies 18 < r' < 22), which weakens the shear 
signal, but this is more than compensated by the extremely large sample size of ~ 16 million 
fgg/bgg pairs. These data indicate that halos of typical galaxies continue an isothermal profile to 
a radius of at least 260 h^ 1 kpc. The conversion from shear amplitude to mass density in the F00 
data is still reliant upon photometric redshift distributions for the fggs and bggs. 
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We present here a detection of galaxy-galaxy lensing around fggs with redshifts determined 
by the Las Campanas Redshift Survey (Shectman et al. 1996) We thus have knowledge of the fgg 
luminosity and distance, and the impact parameter of the lensed photons, in physical units (i.e. 
kpc rather than arcsec), information that is normally lost in projection. The bgg's (R < 23) are 
at magnitudes within the reach of current pencil-beam redshift surveys, and hence their distance 
distribution is well enough determined to introduce negligible uncertainty in the calibration. We 
are therefore in theory able to measure the galaxy luminosity L, halo mass M, and surface mass 
density £ in physical units without ambiguity (save the scale factor Hq). In practice our S/N is too 
low to constrain the extent of the galaxy halo meaningfully, so we combine the F00 profile shape 
with our normalization to obtain a mass profile in absolute units. Combining the resultant M/L 
with the well-determined luminosity function of the LCRS galaxies (Lin et al. 1996) we obtain the 
galactic contribution to fi. We assume an Einstein-deSitter Universe for the values presented here, 
and use Hq = 100/ikms -1 Mpc -1 . Because the fggs are relatively nearby, and the bggs are far 
behind them, the dependence of the results upon the cosmological parameters and A is negligible 
(< 2%). 



2. Observations and Reduction Methods 
2.1. Observations 

We observed 36 square degrees of sky in R band using the Big Throughput Camera (BTC) 
(Wittman et al. 1998) on the CTIO Blanco 4-meter telescope over three runs from December 1996 
through March 1999. Most areas of the target fields are covered by three dithered 420-second 
exposures. Observations were made within the Las Campanas Redshift Survey fields and yield 
« 45,000 well-measured bggs per deg 2 , with typical seeing of l'/l FWHM and RMS sky noise 
of ~ 28i?magarcsec~ 2 on C//43 pixels. Images are debiased and flattened using standard IRAF 5 
routines. Rather than sum images, we analyze each 300s exposure separately and average the 
resulting galaxy measurements. This allows us to adjust for exposure-to-exposure seeing variations 
as well as to avoid or isolate various systematic errors. Observations of i?-band standards from 
Landolt (1992) yield photometric solutions with worst-case RMS errors of 0.05 mag. 



2.2. Shape Measurement 

Galaxy shape measurements are critical to weak lensing. Detailed descriptions of our methods 
are presented in Bernstein et al. (2000) and Smith (2000). Here we summarize the processing steps. 



5 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association 
of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science 
Foundation. 
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Initial detection, photometry, and size estimation is done by SExtractor (Bertin & Arnouts 1996). 
This information is then fed to our shape measurement code elliptomatic, which determines 
cllipticities from Gaussian-weighted second moments of each object. The Gaussian weights are 
elliptical; the weight ellipticity is iterated to match the object's ellipticity, which produces an 
unbiased estimate of the object shape. Stellar objects are selected from the catalog, and a position- 
dependent convolution kernel is created which, when applied to the original image, makes the stellar 
PSF round everywhere in the image (Fischer & Tyson 1997; Bernstein et al. 2000) This corrects 
all objects for shape errors induced by telescope tracking, seeing, or linear CCD charge transfer 
inefficiencies. Optical distortions are significant across the wide field of the BTC. Observations 
of astrometric standard fields yield a distortion map, and object shapes are corrected for this 
distortion analytically, yielding the observed (post-seeing) ellipticity components e° and e%. The 
observed shapes are then corrected for the circularizing effects of the PSF to give an estimate of 
the pre-seeing ellipticity e\ and e l 2 which would be observed in an image with perfect resolution. 
If this correction is larger than a factor of 5, the object is discarded. The measured ellipticities of 
a given object on different exposures are averaged to give the mean pre-seeing ellipticity for that 
object. An estimated uncertainty a e due to image noise is also calculated for each object. 

When considering the distortion induced by a given fgg, the bgg ellipticity is rotated to com- 
ponents e+ and e x , where the former is positive for an image oriented tangentially to the fgg. 

The algorithms used herein to correct for PSF anisotropy and circularization are known to 
leave systematic residuals at the ~ 0.5% level in the shapes of the galaxies. These have little effect 
upon the galaxy-galaxy lensing measurements since one is measuring (e+) in annuli about the fggs, 
and a systematic distortion of constant orientation will cancel to first order when integrated around 
the annulus. We will demonstrate the absence of significant systematic errors below. 



2.3. Determination of Distortion from Ellipticities 

A given bgg with source-plane tangential ellipticity e s + , when subjected to a weak (tangential) 
distortion S, will in the absence of noise be measured to have image-plane ellipticity 

e\ = e% + 8[1 - (e s + ) 2 ]. (1) 

Since the bgg orientations are assumed isotropic, (e+){, = 0, a useful unbiased estimator for the 
distortion is 

8 = (e\) b /K, TZ=(1- {e\) 2 ) b = 1 - <t 2 sn , (2) 

where <jg N = ((e t + ) 2 )i ) is often called the shape noise. The responsivity 7Z is similar to the mean 
"shear polarizability" of Kaiser, Squires, & Broadhurst (1995). In the presence of measurement 
error a e upon each ellipticity, we may wish to apply some weight function w(e,a e ) to each bgg 
shape. The estimator for local distortion becomes (Bernstein et al. 2000) 

i = < 3 > 



-5 - 



E fc ^(l-/4) + ^^(l-/e 2 ) 

n = — l - 



Eb w 



(4) 



The parameter / describes the fraction of the variance in measured e*,_ for the bgg population 
which is attributable to shape noise. Any weight function is permissible if it depends only upon 
the magnitude, not the direction, of the ellipticity e. We use the following weight, which can be 
shown to minimize the variance of the estimator under certain conditions (Bernstein et al. 2000): 



,(e.. e)= 1 -^»: + ^W-'> . (6) 



Alternative choices of weight function lead to results that are indistinguishable within the noise. 

We have tested the accuracy of the distortion-measurement algorithms using simulated data 
frames. An image of artificial spiral and elliptical galaxies is distorted by a circular isothermal 
potential, then subjected to degradation by Gaussian seeing and noise similar to those in our 
images. The artificial images are then run through the software pipeline, and the above equations 
used to extract the distortion 5 as a function of radius from the putative lens. The results are 
plotted in Figure 1, which illustrates that the input distortion is recovered to an accuracy of 5% or 
better as long as the distortion is indeed weak. To insure that the calibration of PSF circularization 
is correct and robust, we repeat the simulation with the angular sizes of the bgg's arbitrarily reduced 
by 30%. The input distortion is again recovered to the same accuracy. 

An additional complication of the real data is that the bgg's are distributed in redshift but 
the distances to individual galaxies are not known. Pencil-beam redshift surveys give the joint 
distribution N(m, z) of the background population. We will assume that the responsivity 1Z and 
weight w of the bgg's are uncorrelated with bgg redshift z b at a given bgg magnitude m b . Then 
the 5 quantity in Equation (3) is a scaled estimator of the distortion 5oo that would be measured 
for a bgg population at infinite redshift: 

<*> = s j^"^'^ (7) 
l^b w 

I^(. D fb/D ob )N(m b ,z b )dz b 

F D (z f ,m h ) = —Z ^oo— — . (8) 

J N(m b ,z b )dz b 

The -D's are angular diameter distances between observer, foreground, and background redshifts. 
The quantity 5oq is thus measurable as long as the distribution of the background weight vs redshift 
is known. The redshift distribution of our bgg sample is sufficiently well known from deep redshift 
surveys that the uncertainty in the factor Fp does not dominate the error budget, as will be 
discussed in detail below. 
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Distortion Calibration Simulations 
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Fig. 1. — Results of shear calibration with simulated data: a simulated isothermal distortion is applied to 
artificial galaxy images, which are then measured via the standard measurement pipeline. The ratio of the 
measured to the applied distortion is shown for our standard weighting scheme (solid line) and two alternative 
weighting schemes (dashed lines). Ideally the ratio is unity to order e, i.e. will fall between the two dashed 
curves. The measured values are within ~ 5% of the applied distortion in the limit of weak lensing, and are 
nearly independent of weighting scheme. The results are not visibly changed when the artificial galaxy sizes 
are reduced by 30%. 
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2.4. Mass Density and Scaled Distortion 

The physical quantity of interest is the mean azimuthally averaged surface mass density £ l (R) 
at radius R about galaxies of luminosity L. The measured quantity is the lensing distortion Soo 
induced upon a distant bgg z^, which is related to the surface density by 

5M L, z f ) = [E L (< R) - Z L {R)] ^D of (9) 

Here Xl(< R) is the average surface mass density interior to R. We will wish to combine data 
from fggs which span a range of Zf, and we additionally will bin data over finite ranges in impact 
parameter R, and fgg luminosity L. To compensate for the spread in these quantities, we express 
the distortion S OQ (R, L, zf) for a given fg/bg pair in terms of a scaled distortion 

<5* = S(R = Ro,L = L*,zf = z*,z b = oo) (10) 
= S 00 (R,L,z f )/[F L (L)F R (R)F z (z f )}. (11) 

Our canonical luminosity is the in the Schechter-function parameterization of the LCRS galaxy 
luminosity function, which is absolute i?-band magnitude M* = —20.3 + 5 log ft,, or = 8Ah~ 2 x 
1O 9 L (Lin et al. 1996) z* is taken to be 0.1, the median redshift of the LCRS sample; and Rq is 
taken as lOO/i -1 kpc if not otherwise specified as the center of some radial bin. The scaling factor 
for redshift is known: 

F z (z f ) = D of (z f )/D of (z*) (12) 

Scaling the measured shear to a standard luminosity and radius requires some model for how 
the surface density profile varies with these quantities. With a sufficiently large galaxy sample, 
one could simply measure the function £l(< R) — ^l(R) in narrow bins of L and R, and non- 
parametrically map this density contrast as a function of these two variables. Our samples are, 
however, small enough that our bins must cover a substantial range of L and R, so we will adopt 
the assumption that the surface density is described by an isothermal profile with a circular velocity 
v c that is a power-law function of luminosity: 

^l(R) = S L (< R)/2 = v c (L) = Vif {L/L^I\ (13) 
The scaling functions are then very simple: 

F L = {L/L*f, F R = R /R. (14) 

The canonical distortion is 

= 2n4^4^ = 0.0071 ( "* , ) 2 ( ^ ) " (15) 
c 2 Rq V 200 km s- 1 / \100h~ 1 kpc J V ; 

The high-S/N data of F00 produce a mean distortion profile that is fully consistent with the 
isothermal profile at the radii R < 200h~ 1 kpc in which we measure S. The power-law dependence 
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of £ upon galaxy luminosity L is suggested by the Tully-Fisher and Faber-Jackson relations, which 
lead us to choose the value (5 = 0.5. These relations have, however, been tested only at radii of 
~ 30/i _1 kpc or smaller — the F00 study cannot test this since the fgg luminosities are not known 
individually. Indeed Zaritsky et al. (1997) were unable to detect any dependence of halo mass upon 
luminosity in the spiral satellite study. We will also calculate 5* under the assumption (3 = 1, the 
mass-traces-light value. In §3.5 we will use our own data to bound (3. 

The goal of this paper is an absolute determination of 5*. Our estimator for 5* combines 
Equations (3), (7), and (10). In addition the variance of the estimator is minimized by adding a 
factor of {FlFrF z Fd) 2 to the weight function. We thus obtain 

~ _ 1 E ftb wF L F R F z F D e% 

nZ f , b w[F L F R F z F D ]i [Lb) 

Standard propagation of errors, along with the substitution Var(e + ) = (e+), yields the variance of 
the above estimate: 



3. Results 

Our foreground sample consists of all 790 galaxies from the LCRS redshift survey with 0.05 < 
Zf < 0.167 that lie within our imaged areas. Background galaxies are required to be well-measured 
on at least two exposures, have magnitude 20 < R < 23, and correction for PSF circularization at 
most a factor 5 (this is effectively a lower limit on angular size). This leaves « 450,000 fgg/bgg 
pairs with impact parameters of 15/i" 1 kpc < R < 480/i _1 kpc. 



3.1. Background Galaxy Redshift Distribution 

An absolute determination of mass densities requires an estimate of the distance factor Fr>(zf, mi,) 
from the weighted redshift distribution of the background galaxies. The Caltech Redshift Survey of 
the Hubble Deep Field and its flanking fields (Cohen et al. 2000, CRS) provides a virtually complete 
redshift survey to R < 23. We combine the CRS N(m, z) data into 0.5-mag bins to construct the 
function Fd- This factor is above 0.6 for nearly all the fgg/bgg pairs in our data, and is by def- 
inition < 1. This limited range of variation is a consequence of the substantial magnitude (hence 
distance) gap between our foreground and background samples, and makes our calibration very 
stable against perturbations in the assumed N(m, z). Thus while the HDF field may be deficient 
in nearby galaxies, we do not expect a significant effect upon our results. 

To estimate the uncertainty in <5* due to Fd errors we turn to the Canada-France Redshift 
Survery (Crampton et al. 1995, CFRS), which is nearly complete for a sample of « 600 galaxies in 
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5 fields defined by Iab < 22.5. We have imaged two of these fields in R band using the 2.4-meter 
Hiltner telescope at MDM Observatory to obtain an N(m,R,z) distribution of the CFRS sample. 
We can either use the measured iV(m#, z) in these two fields, or use them to construct a mean 
(R — Iab) vs (Vab — Iab) color relation for the galaxies, and apply this relation to the full CFRS 
catalog to mimic an .R-band observation. The values of 5* determined by the two methods differ 
by only 1%. 

Using the measured CFRS N(m, z) changes 5* by ±5% or less relative to the use of the CRS. 
The CFRS is substantially incomplete in two senses: first, the Iab < 22.5 cutoff of the CFRS 
catalog misses an unknown number of the bluer R < 23 galaxies. These bluer galaxies are likely 
nearer, on average, then the measured R ss 23 sample. Second, the CFRS redshift measurements 
are incomplete even for the Iab < 22.5 sample; this incompleteness is as large as 30% in the 
22.5 < R < 23.0 bin. These galaxies are likely to be biased toward the 1 < z < 2 range where 
redshift measurement is difficult. If we place all the unmeasured galaxies at z = 2, 5* is still only 
6% lower than the CRS result — and it is likely that inclusion of the missing blue R 23 galaxies 
would push this number back toward the CRS result. 

We conclude that calibration errors due to uncertainties in the bgg redshifts are < 5%. 



We plot the scaled distortion <5* as a function of proper distance from the lens galaxy center in 
Figure 2. The tangential alignment is clearly detected to ~ 150/i _1 kpc. As a test for systematic 
contamination, we calculate the azimuthal average of e x - The e x averages are consistent with zero 
distortion, as they must be if the distortion is due to weak lensing (Stebbins, McKay, & Frieman 
1996), giving us confidence that systematic shape distortions have minimal effect upon our <5* 
determination. 

Our data are consistent with an isothermal profile but are too noisy to provide any useful 
constraint on the profile shape. No useful lower limit, for example, is implied upon a possible outer 
truncation radius R max for an isothermal halo. The high-S/N data of F00, however, do provide 
useful constraints to the truncation radius which we will use below. 

The most accurate determination of 5* is made by fitting all the data with 15/i _1 kpc < R < 
240/i~ 1 kpc to the isothermal model. For the two trial values of (3 we obtain 



3.2. Measured Halo Properties 



<5* = 



{ 



0.0056 ± 0.0014 = 0.5 
0.0048 ± 0.0010 = 1.0 



(18) 



In either redshift model the lensing signal is detected at signal-to-noise ratio of 4-5. For comparison, 
the <5 X values computed over this broad bin are (1. ± 11) x 10~ 4 ((3 = 1), consistent with zero within 
the calculated uncertainties. 
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Fig. 2. — The square symbols show the measured distortion 6* (scaled to an L* galaxy at z/ — 0.1 and 
Zb = oo ) vs impact parameter. The dashed line is the best-fit isothermal model. The cross symbols (x) 
show the binned means of the S x component of shear, which should (and does) vanish if the distortions are 
due entirely to lensing. This plot depicts results assuming M oc L* 3 with (3 = 1, and the total significance of 
the detection is 4.5a. The results for (3 = 0.5 are similar. 
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3.3. Sources of Error 

Possible sources of error in our determination of <5* include: 

1. Noise: the largest source of error is the random shape noise and measurement noise in the 
ellipticities of the bggs. This contributes an uncertainty of 20-25% (la) to <5*, and is the error 
quoted in Equation (18). 

2. Background Redshifts: As discussed above, the R < 23 redshift distribution is known well 
enough to reduce calibration uncertainty to ±5% or less. 

3. Distortion Calibration: The accuracy with which our software measures the galaxy shapes, 
corrects to the pre-seeing shapes, and calculates a distortion has been tested by the simula- 
tions. This accuracy appears to be ss 5% or better. 

4. Mis-measurement: One could imagine that the shape measurements of background galaxies 
could be biased by the wings of the foregound galaxies. At radii beyond 30/i _1 kpc, however, 
the light from the foreground galaxies is far below the noise of the images. 

5. Luminosity Scaling: The measured (5* depends to some degree on the index (3 scaling lumi- 
nosity to mass [Equation (13)]. We will attempt to measure (3 below. 

6. Satellite Galaxies: Our sample of "background galaxies" will be diluted by faint galaxies 
that are physically associated with the foreground galaxies. (5* should be adjusted upward 
by a factor equal to the fraction of bgg's that are in fact satellite galaxies. F00 measure 
this contamination to be about 15% at a projected radius of ~ GO/i^ 1 kpc, where our data 
are centered. Our fgg sample is similar to that of F00, but the bgg magnitude ranges differ: 
18 < r 1 < 22 for F00, and 20 < R < 23 for us. Our fainter bgg sample should have a 
substantially smaller fraction of physically associated galaxies since background counts rise 
much more quickly (factor of 2.5 per magnitude in R, (Tyson 1988)) than the dwarf galaxy 
luminosity function (factor of 1-1.5 per mag). We detect this effect in our own data, but it 
is difficult to measure precisely. We estimate, however, that our contamination factor should 
be below the 10% level. 

7. Projected Neighbors: The distortion S(R) measures the surface mass density (above the cosmic 
mean density) projected within radius R of a fgg. As galaxies are clustered in space, a circle 
of radius R about a fgg actually includes mass from the halos of excess neighboring galaxies. 
The strength of this effect can be estimated by convolving the galaxy-galaxy correlation 
function £(r) f» (r/8h~ 1 Mpc)~ 1 - 8 with an isothermal mass distribution. At our radii of 
R ~ 100/i -1 kpc the effect is quite small (« 4%) and can be ignored. Beyond 200/i~ 1 kpc 
the effect rises above 10% and must be taken into account (see F00). Note that this effect 
opposes that of satellite "dilution" mentioned previously, but both tend to flatten the run of 
5 with radius. 
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In summary, the errors in 5* are dominated by random noise at ±20%. There are 4 calibration 
uncertainties each at < 5% level, which are known to work in dissimilar directions, so the overall 
calibration uncertainty should be below 10% (save the uncertainty in (3), and will be ignored. 



3.4. Total Mass in Galaxy Halos 

The halo mass profiles are consistent with an isothermal profile to the ~ 200h~ 1 kpc radii we 
probe; the total mass associated with each galaxy depends upon the extent of these isothermal halos. 
If the surface density of a galaxy follows the isothermal profile in Equation (13) to a projected radius 
Rmux, beyond which it is abruptly truncated, then the total mass of a typical galaxy at luminosity 
L is 

M(L) = M*(L/L*f, (19) 
M * = ^ v * Riaax ( 2 °) 
- 2 .06xl0»™,y^-). (2!) 

Our data do not constrain i? max , so we turn to the SDSS data. The F00 distortions are fit to a 
model in which all galaxies are truncated at some angular size s; F00 obtain a 95% CL lower limit 
of s > 140" for their sample. 6 The estimated redshift distribution of fggs in the F00 sample has a 
weighted mean angular diameter distance of 0.125c/ Hq, so a simple estimate of the physical scale 
of the truncation radius is s > 260/i _1 kpc at 95% CL. The F00 model ignores possible variations 
of halo extent with galaxy luminosity, but such considerations are not likely to influence our result 
because the mass integrals and the F00 data are both dominated by galaxies near L*. We will 
therefore adopt from F00 the result 

-Rmax > 260/j- 1 kpc (95% CL). (22) 

Taken with Equations (18) and (21), this implies that the total mass and mass-to-light ratio 
for an L* galaxy contained within projected radius of 260/i~ 1 kpc are: 

M{<n)IL . {^tZZfC 5:13 

Note that the LCRS luminosities are isophotal, and Lin et al. (1996) recommend raising the lumi- 
nosity estimates (hence lowering M/L) to account for flux excluded from the apertures. There has 
also been no account taken of evolution since z ~ 0.1, though a K-correction is made. 



6 The fitting function in F00 is not quite a cylindrical truncation at s, but the integrated mass is the same. 
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Using the Schechter-function parametric fit to the LCRS galaxy luminosity function from (Lin 
et al. 1996) and the mass formula in Equation (19), the mean mass density in galaxy halos is 
M*0*r(l + a + 0), where <p* = 0.019/t 3 Mpc -3 and a = -0.70 are the fitted Schechter-function 
parameters. As a fraction of the critical density, the density in LCRS halos is 



These constitute a lower limit to f2 since the matter associated with galaxies may well extend 
beyond 260/i _1 kpc. Furthermore we have only inventoried matter which is associated with galaxies 
included in the LCRS sample. Low-surface-brightness galaxies excluded from the LCRS would 
increase Q to some extent, as of course would any form of matter which is not spatially correlated 
with the centers of LCRS galaxies. 



The largest uncertainty in our calculation of f2 is the exponent (3 describing dependence of 
galaxy mass upon luminosity. This parameter is normally assumed to be near 0.5 based upon 
dynamical measures in the luminous parts of galaxies, but it is purely an ansatz for mass at radii 
beyond 30/j _1 kpc. We have divided our foreground sample into 3 bins of luminosity and fit an 
isothermal distortion profile to an annulus 15/i -1 kpc < R < 240ft, -1 kpc for each bin. In Figure 3 
the dependence of distortion (mass) upon L is readily apparent. This is in contrast to the Zaritsky 
et al. (1997) satellite galaxy study, which did not see a significant correlation between luminosity 
and halo mass at 200/i _1 kpc. The galaxies certainly "know" about the density of the halos in 
which they are embedded. 

The dependence of 5 upon L is readily fit by the Equations (13) for the mass-traces-light 
value of (3 = 1, (x 2 /v = 1-16/2). The Faber-Jackson value of (3 = 0.5 is somewhat disfavored 
(x 2 l v = 4.93/2, Q = 0.09). Treating (3 as a free parameter in a fit to the three bins yields a 95% 
CL range of 0.6 < (3 < 2.4. This result is intriguing and demonstrates the power of weak lensing 
to elucidate the relation between luminous and dark matter. The traditional (3 = 0.5 model is not 
strongly excluded, however, so we will await further data before drawing conclusions. 

A similar marginal difference is seen between the masses of absorption-spectrum and emission- 
spectrum galaxies in the foreground sample, but meaningful results must await better statistics. 




> 0.23 ±0.06 (0 = 0.5) 

> 0.16 ±0.03 ((3=1.0) 



(25) 



3.5. Mass vs Luminosity 
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Fig. 3. — The measured distortion 5 (at R = 60/i _1 kpc) is plotted for three bins of foreground galaxy 
luminosity L. Dependence of distortion (and hence mass) upon L is clearly detected. The dashed line 
is the best-fit mass-traces- light model (M oc L), which fits the data well; the dotted line is the best-fit 
Faber- Jackson model (M oc L ' 5 ), which is only marginally acceptable. 
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4. Comparison to Other Determinations 

4.1. Comparison with Other Galaxy-Galaxy Lensing Estimates 

We place our result in context first with previous measures of galaxy-galaxy lensing. Most 
previous works were parameterized by the halo circular velocity of an galaxy in the isothermal 
portion of the halo, which we have called u*. Our measure of 5* gives 1-sigma limits on v* of 

f 176 ±22 km s" 1 (0 = 0.5) 
v* = \ i 26) 

\l64±20kms _1 (0 = 1.0). 

This is in reasonable agreement with the determinations of Brainerd, Blandford, & Smail (1996) 
(v, = 220±80kms~ 1 , 90% CL), Hudson et al. (1998) (v, = 210±40km S - 1 ) and of F00 (200-280 
km s -1 , 95% CL), given that the disparate definitions of L* have been used (the LCRS value is on the 
low side). We emphasize, however, that ours is the first determination of this quantity to actually 
measure the luminosities of the foreground galaxies to give a result which is not heavily dependent 
upon assumptions about the luminosity and redshift distributions of the galaxy populations. 



4.2. Comparison with Other Dynamical Estimates 

Zaritsky & White (1994) use the satellite galaxy data to constrain the mass M200 enclosed 
within a 200/i _1 kpc sphere of their typical spiral to be 1.1— 2. Ox 10 12 /t _1 M@. Our results imply a 
mass within 200/j- 1 kpc of (1.25 ± 0.25)(L/L*) x lO 12 /^ 1 M Q for = 1, or (1.44 ± 0.36) y/L/L* x 
10 12 /i _1 M Q for (3 = 0.5. There is general agreement, but a more specific comparison is not 
possible because Zaritsky k, White find no dependence of mass upon primary luminosity (0 « in 
our terminology) while we find a relatively strong dependence ((3 > 0.5). 



4.3. Comparison to the Tully-Fisher Relation 

It is of interest to compare the circular velocity of the halo to the circular velocity of the gaseous 
disk. Tully et al. (1998) present a calibration of the Tully-Fisher relation (TFR) for Ursa Major and 
Pisces cluster spirals in the R band. The slope of their .R-band TFR is equivalent to (3 = 0.62. But 
they also find that the extinction is stronger in brighter galaxies, and the LCRS magnitudes are not 
corrected for extinction, so the mass-luminosity relation in the LCRS data would appear somewhat 
shallower, (3 = 0.58, for a typical galaxy at ~ 60° inclination. For this value of (3, the lensing data 
suggest a circular velocity v* = 176 ± 22kms _1 at an effective radius of ss 60ft, -1 kpc. The TFR 
calibration yields a disk circular velocity of 160 km s -1 for a typically-inclined spiral at the LCRS 
value of L*. This consistency suggests that the isothermal mass profile observed at ~ 10/i -1 kpc 
in the disk matches fairly well onto the halo profile at ~ 100/i _1 kpc. Even if galaxies have perfect 
isothermal profiles, we might expect the lensing v* to be slightly higher than the TFR value (for 
late-type spirals) because of the inclusion of early types in the lensing foreground sample. 
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5. Conclusions 

The relations between galaxy luminosity and halo mass determined herein provide basic con- 
straints to models of galaxy formation and evolution. We find, perhaps surprisingly, that the 
halo circular velocity at 60-100/i _1 kpc radii are similar to the disk circular velocities measured at 
< lO/i^ 1 kpc, extending the "disk-halo conspiracy" to larger radii than previously known. There 
are hints in the data, however, that this conspiracy breaks down in the sense that halo masses vary 
more sharply with galaxy luminosity then would be indicated by the Tully-Fisher or Faber-Jackson 
relations. 

The matter density £l m has been measured by many methods in the past, and there is a 
currently fashionable consensus that Q m pa 0.3. Our results (Equation [25]) are in agreement with 
this consensus given that our method produces only a lower limit. Our limits are, however, close 
to the consensus value, especially if (5 pa 0.5. Our lower limit on Q, is determined by a direct 
inventory of the matter in the Universe; the only untested assumption is that the formula for 
gravitational deflection of light can be extrapolated from Solar-System scale (where it has been 
verified) to galactic scales. Other determinations of £1 rely upon assumptions about the nature 
of our Universe. Agreement between our direct inventory and these other measures serves as 
valuable and unique verification that these assumptions — often fundamental to our current view of 
cosmology — are correct. 

The more precise estimates of £l m include the following: Dynamical studies of galaxy clusters, 
combined with the assumption that the M/L ratio for clusters is universal, yield f2 m = 0.24+° og 
(Carlberg et al. 1996). Another cluster-based means to f2 is to measure the fraction /& of the cluster 
matter that is baryonic, assume that the cluster ft> is universal, and then use the baryon density 
£lb from Big Bang nucleosynthesis to estimate Q, m = Q^/fb- Mohr, Mathiesen, & Evrard (1999) 
conclude that Q m < 0.30 ± 0.04 (for h = 0.7 ± 0.1). 

Other estimates of measure its effects upon global geometry or the growth rate of fluctuations 
in the Universe. In the former category, the combination of high-redshift supernovae measurements 
and the location of the cosmic background Doppler peaks measured by MAXIMA- 1 constrain 
0.25 < Q m < 0.50 at 95% CL (Balbi et al. 2000). Peculiar velocity studies determine the parameter 
Pi pa Sl 6 /t> to be pa 0.5 ± 0.05 (Strauss & Willick 1998), with b a galaxy bias parameter. Taking 
6 pa 1 gives a value of £1 pa 0.3, again slightly above our halo density. Conversely our lower limit on 
f2 implies a bias parameter for IRAS galaxies b > 0.7. 

Thus our results may be taken as independent evidence that these underlying assumptions are 
correct: cluster M/L and baryon fractions can be taken as representative; general relativity does 
properly describe the geometry of the Universe to z pa 1000 given our directly measured matter 
content; and large-scale motions are consistent with the predictions of gravitational instability. 

Conversely, we may take the £l m pa 0.3 value as truth and conclude from these lensing results 
that most of the matter in the Universe is indeed contained within 260/i _1 kpc of normal galaxies. 
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This implies that galaxies missed by the LCRS survey — e.g. low-surface-brightness galaxies or 
exotic dark galaxies — cannot be the dominant reservoirs of mass in the Universe. Likewise there 
cannot be a dominant matter component which is not clustered around galaxies. Most of the 
mass in the Universe appears to have been located, though of course its material nature remains 
unknown. 

The accuracy of our results is currently limited by the number of foreground/background pairs 
in our survey. This will be remedied in dramatic fashion by the full Sloan Digital Sky Survey, which 
will obtain imaging to r' < 22 around roughly three orders of magnitude more foreground galaxies 
than used in this study with a concomitant increase in the number of pairs. This will allow for 
non-parametric mapping of the shear profile as a function of radius, galaxy luminosity, and galaxy 
type, with much lower noise than in our LCRS lensing survey. This weak lensing information will 
be unique in unveiling the relation between the various classes of galaxies and the dark halos in 
which they reside. With the random errors reduced, closer attention will have to be paid to some 
of the < 5% corrections we have ignored (as in F00), but these are tractable. 

Our estimates of M* and are lower bounds since the profile of the halo has not yet been 
determined beyond 260/i _1 kpc. Beyond this radius it becomes inappropriate to speak of "the halo" 
of a galaxy, as the mass profile of a given galaxy becomes inseparable from the mass associated 
with neighboring galaxies. More stringent limits on 0, will thus require a more elaborate formalism 
for describing the effects of galaxy correlations upon the weak lensing measurements. 

This work was supported by grant AST-9624592 from the National Science Foundation. We 
gratefully acknowledge the assistance of J. A. Tyson and David Wittman with the construction 
and operation of the BTC camera and their contributions to the data pipeline; Tim McKay and 
Erin Sheldon for continued discussions of SDSS lensing results; and the staff of CTIO for their 
enthusiastic support of the BTC observations. 
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